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Параметрическое самовозбуждение динамической системы резания” 


В. Л. Заковоротный, Фам Тху Хыонг 
(Донской государственный технический университет) 


Рассматривается проблема обеспечения устойчивости траекторий формообразующих движений инструмента 
относительно заготовки с учётом параметрического самовозбуждения. Анализируются факторы, вызываю- 
щие периодические изменения параметров в динамической системе резания. Они связаны с пространствен- 
ной анизотропией свойств упругости в подсистеме обрабатываемой заготовки, с вариациями припуска по 
периметру вращения заготовки, с кинематическими возмущениями со стороны механической части приводов 
исполнительных элементов станка, с периодическими процессами в зоне резания. Приведена обобщённая 
динамическая модель системы с периодически изменяющимися параметрами. Выполнено изучение влияния 
периодически изменяющихся параметров на устойчивость траекторий. Раскрыты особенности потери устой- 
чивости, специфические для динамических систем резания. В частности, показано, что за счёт параметриче- 
ских эффектов по мере увеличения скорости резания при возрастании частоты вращения шпинделя всегда 
существует критическая частота, при которой система теряет устойчивость. 
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Введение. Проблема обеспечения устойчивости траекторий формообразующих движений ин- 
струмента относительно заготовки является одной из самых важных при выборе технологических 
параметров и конструктивных особенностей подсистем инструмента и заготовки, взаимодейству- 
ющих через процесс резания [1—3]. Существующие представления о потере устойчивости бази- 
руются на использовании в качестве математических моделей систем линейных или нелинейных 
дифференциальных уравнений динамики. Однако практика показывает, что существующие пред- 
ставления противоречат многим известным экспериментальным данным. Например, исходя из су- 
ществующих представлений, по мере увеличения скорости резания запас устойчивости в системе 
должен возрастать. Однако практика показывает, что это справедливо при относительно малых 
частотах вращения шпинделя. При увеличении частоты до некоторой пороговой величины всегда 
наблюдается увеличение вибраций, свидетельствующих о потере устойчивости траекторий. В 
настоящей статье, основываясь на известных представлениях о динамической системе резания, 
рассматривается практически не анализируемый ранее вопрос о самовозбуждении системы за 
счёт периодического изменения параметров в математическом описании динамики системы. 
Оснований для учёта переменности параметров в динамической системе резания множе- 
ство. Во-первых, свойства упругости в подсистеме обрабатываемой детали не являются симмет- 
ричными. Это связано, по крайней мере, с изменениями радиальной жёсткости детали, закреп- 
лённой, например, в трёхкулачковом патроне. Во-вторых, при обработке заготовки, имеющей пе- 
риодические изменения припуска, параметры динамической характеристики процесса резания 
также являются периодически изменяющимися. В-третьих, траектории исполнительных элементов 
станка всегда являются периодически изменяющимися за счёт кинематических возмущений, зави- 
сящих от точности и конструктивных несовершенств механической части приводов. В-четвёртых, 
при обработке большинства материалов имеет место периодическое изменение сопротивления 
резанию, связанное с периодическим формированием поверхности скольжения в области первич- 
ной пластической деформации. Наконец, многие технологические режимы принципиально имеют 
периодические изменения параметров, например, процесс фрезерования. Все перечисленные 
факторы вызывают, по крайней мере, периодические изменения суммарной жёсткости в уравне- 


* Работа выполнена в рамках инициативной НИР. 
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ниях взаимодействующих подсистем. В таких системах, как известно из теории колебаний, воз- 
можно параметрическое самовозбуждение системы, зависящее от частоты периодического изме- 
нения параметров. 

Обобщённая математическая модель динамической системы резания как системы с 
периодическими изменениями параметров. Для изучения условий самовозбуждения дина- 
мической системы резания, параметры которой являются периодически изменяющимися функци- 
ями времени, необходимо рассмотреть обобщённую модель системы. Такая модель должна учи- 
тывать главные свойства системы, но не быть излишне сложной. Ограничимся случаем, когда 
рассматриваются периодические изменения элементов матрицы жёсткости. Кроме этого, примем 
во внимание гипотезы, используемые ранее при анализе динамики системы с постоянными пара- 
метрами [2]. Они включают три положения: силы, формируемые в зоне резания, пропорциональ- 
ны площади срезаемого слоя; они имеют неизменную ориентацию в пространстве, задаваемую 
угловыми коэффициентами х,, / =1,2,3; модуль силы запаздывает по отношению к деформаци- 


онными смещениями инструмента относительно заготовки. Таким требованиям соответствует сле- 
дующая модель, рассматривающая динамику системы в вариациях относительно точки равновесия: 








Фу ау ау ау | 
т се +1, > +В, Е +В. а. + [с,1 + А, (9,6) у, + 


+ [2 + А, › п (9,Ё + 6, ) | уч [2 ++ рхиЁ. (1+ нп (0,2) | У; +РХ,5. 02, =0; 








ау ау ау ау 
т 2 +1, > Е +1, › ОЕ +2, > — + [Се + А, > 5т (0.6 +6, , ) у, + 


+] С,> + А, яп(® +6, у» + [@ + Хх, „о (1+5 (9,2) Уз+РХ,5, 02, =0; (1) 
а?у ду. ау ау 
а ЗЫ 
+ - ++ рХ.Ё,, (1+5 (9,2) У; + РХ.5ро2, = 0; 
[@74 
Т, Е +2, =У,, 


где у={у,,У,, у — вектор упругих деформационных смещений инструмента в простран- 
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стве, мм; р — давление стружки на переднюю поверхность инструмента, параметр, завися- 
щий от физико-механических свойств обрабатываемого материала, геометрии инструмента и 


технологических режимов, кГ/мм?; $0’ во — Соответственно величины подачи и глубины 


ро 
резания, не учитывающие упругие деформационные смещения; т= [1 5,К =1,2,3; 
т, к =0, при: $ #К, й= [1 ‹ 15, К =1,2,3, с= Е 5,К =1,2,3 — соответственно матри- 
цы инерционных, диссипативных и упругих коэффициентов (они являются симметричными и 
положительно определёнными); 7. — постоянная времени, моделирующая запаздывание ва- 


риаций модуля силы резания по отношению к деформационному смещению у/, с. 
Таким образом, при установившемся состоянии в точке равновесия 2, = у, =0. Кроме это- 


го, в (1) учитывается пространственная анизотропия упругости системы в виде матрицы: 
А. 5т(9#) А 9т(0Ё+0,,) 0 
А=| А, эп (96 +60,,) А, эт(9Ё+60,,) 0). (2) 
0 0 0 
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В (2) отсутствует переменная составляющая по направлению у. , так как кинематические возму- 


щения по направлению оси вращения шпинделя формируются на основе следующего интеграль- 


ного преобразования: 
Е 


у = Ги (Е) 4%, 


ЕТ 


где и, (Е) — переменная составляющая скорости продольной подачи, мм/с. 
Нетрудно показать, что в у,(Ё) все переменные составляющие, равные или кратные пе- 


риоду вращения шпинделя Т’, равны нулю. Причём в (1) учитываются переменные составляю- 
щие, обусловленные вариациями координаты суппорта, в направлении, ортогональном к направ- 
лению его движения, а также вариации припуска, периодические с частотой вращения шпинделя 


0 0 рх.Ё,, (1+нэт(9,2)) 
с» (Е) =|0 0 рх,Ё,„, (1+няп(9,2)) |. (3) 
0 0 рх,Ё, , (1+н5т(9,2)) 
В (1) частоты ©, и ©, равны или кратны частоте вращения шпинделя. Поэтому для анализа 


устойчивости системы (1) можно использовать теорию Флоке [4]. При изучении областей устой- 
чивости использовался метод прямого цифрового интегрирования систем дифференциальных 
уравнений, основанный на методе Рунге — Кутты четвёртого порядка. Потеря устойчивости фик- 
сировалась по траекториям решений дифференциальных уравнений. 
Особенности параметрического самовозбуждения в динамической системе резания. 
Система (1) отличается от традиционно анализируемой системы следующими основными особен- 
НОСтТями: 
— в реальной динамической системе резания всегда имеет место изменение параметров, обуслов- 
ленное различными причинами, поэтому, как правило, имеется несколько частотных составляю- 
щих изменения параметров, кратных частоте вращения шпинделя; 
— вариации сил запаздывают по отношению к изменениям упругих деформационных смещений, 
причём запаздывание зависит как от технологических режимов, так и от геометрии инструмента и 
физико-механических свойств обрабатываемого материала. Поэтому даже без параметрических 
эффектов система может потерять устойчивость. Периодические же изменения параметров, как 
правило, дополнительно способствуют потере устойчивости; 
— упругие свойства подсистемы заготовки без процесса резания являются линейными и матрицы 
упругости являются симметричными и положительно определёнными. Поэтому в пространстве 
для моделирования упругости можно на основе матрицы поворота выделить оси коллинеарного и 
ортогонального направлений, которые вращаются при повороте шпинделя при заданной анизо- 
тропии свойств упругости. Это обстоятельство может вызывать параметрическое самовозбужде- 
ние системы, которое дополняет условия самовозбуждения за счёт влияния циркуляционных сил, 
формируемых динамической связью процесса обработки. 

Для дальнейшего анализа приведём некоторые характерные примеры параметрического 
самовозбуждения системы. Для этого рассмотрим различные диаграммы параметрического само- 
возбуждения в плоскости 


П=6/м, |, 
где и, — собственная частота основного колебательного контура системы, Гц; и — частота па- 
раметрического возбуждения, Гц; и — уровень параметрического возбуждения (и<1). 
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Для более наглядного выяснения влияния различных особенностей системы на потерю 
устойчивости, сначала рассмотрим скалярный аналог системы (1). Проанализируем для этого слу- 
чая влияние двухчастотного параметрического возбуждения для системы: 





2 
т@ . +10 сани по +, япза у = —с,2; 
@ (4) 
т Е 
[8 ь 


где т, В, с — соответственно масса, коэффициент демпфирования и жёсткости подсистемы; 
с› — жёсткость процесса резания, кГ/мм ; Г — постоянная времени стружкообразования, с. 
Рассмотрим изменение области параметрического возбуждения системы на плоскости па- 
раметров | п) при варьировании |, от нуля до 0,8 и при заданных значениях параметров ди- 
намической системы с, =500 кГ/мм, Т =0,1мс, №=0,10 кГ/мм. Соответствующие диаграммы 
ране возбуждения системы приведены на рис. 1. 
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Рис. 1. Пример изменения области параметрического возбуждения системы на плоскости параметров (ц,, п) 








для случаев р, =0,0 (левая диаграмма) и р, =0,8 (правая диаграмма) 


Приведённый пример вычислен для системы, имеющей собственную частоту консерватив- 
ной системы @, =150,0 с". Если рассматривать частоту вращения шпинделя равную 


3000 об/мин, то есть 50 Гц, то вторая диаграмма уже находится в области параметрического са- 
мовозбуждения системы. Однако, если рассматривается скалярный случай, то параметрическое 
самовозбуждение системы наблюдается при высоких частотах вращения шпинделя, характерных 
лишь для скоростного резания. Ситуация принципиально меняется, если анализировать вектор- 
ный случай пространственных колебательных смещений. 

Приведём пример условий параметрического самовозбуждения системы для случая, когда 
деформационные смещения инструмента происходят в плоскости. Линеаризованное уравнение в 
вариациях для этого случая имеет вид: 














тб Е с. [1 ия (ше +6, у, +5 -изт( +6, у, ео 
тб Це с -т(шё+6,,)-у,+с-[+нят(шё +6] у, =-Х.с»2., Е 
аа, 
тан, =у, 
где с, =р$»›. 
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В системе (5) дополнительно учитывается, что запаздывание сил по двум ортогональным 
направлениям отличается. Обратим внимание на то, что в системе (5) имеет место вращение угла 
ориентации эллипса жёсткости. Этот факт вызывает дополнительные условия параметрического 
самовозбуждения системы резания. 

Рассмотрим конкретный пример (рис. 2). Исходные данные: матрица инерционных коэф- 


фициентов т =5,63.10* кГ?/мм; матрица диссипативных коэффициентов р = 0,1 кГ/мм; матри- 
ца жёсткости с = 2000 кГ/мм; с, = 500 кГ/мм; первая собственная частота системы @, =150 Гц; 


вектор коэффициентов ориентации сил Х = {Х,, х. = {0,5; 0,8} . Диаграмма построена по отно- 


шению к первой собственной частоте системы без резания. Угол ориентации эллипса жёсткости 
рассматривается при и=0. 


























п п 

3 | | р р 3 Т р Е | 
2.5 р ЯНИС _ И НЯ Еее 2.5 а ис РАЯ о Е у 
НЕ ЕРЕЕЕЕ ЗЕЕ | 
а ны ес ее авы о ыы 
А И: РЕ В ыы 

р. т 
Пра А. и а в. а аа 

Н . : : |] 1 } | | 1 Н 

0 й Е | | 0 | Е Е | 

0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 
а) 6) 


Рис. 2. Изменение области параметрического возбуждения в плоскости (п, и) при варьировании угла ориентации эллипса 


жёсткости: а — а = 0°; б6— а = 45° 


Здесь большое значение в механизме самовозбуждения имеют условия обмена силовыми 
потоками между подсистемами, рассматривающими колебательные смещения в двух ортогональ- 
ных направлениях в плоскости, нормальной к оси вращения шпинделя. Приведённые данные по- 
казывают, что в рассматриваемой системе параметрическое самовозбуждение наблюдается уже 
при частотах вращения шпинделя, равных (400,0—500,0) об/мин. Более скрупулёзный анализ по- 
казывает, что на условия самовозбуждения в этом случае оказывают влияние следующие факторы: 

1) Известно [2], что за счёт реакции со стороны процесса обработки суммарная матрица 
упругости становится несимметричной. Поэтому без параметрических эффектов (р =0) в системе 


естественным образом формируются циркуляционные силы. В зависимости от направления дей- 
ствия циркуляционных сил (сил, ортогональных к направлению деформационных смещений) и 
направления вращения эллипса жёсткости могут формироваться условия, способствующие само- 
возбуждению, или направленные на стабилизацию равновесия. 

2) Рассматриваемая система имеет две резонансные частоты, которые разделены в ча- 
стотной области антирезонансом. При сближении этих частот в системе формируются формы ко- 
лебаний типа биений с частотой, равной разности резонансных частот. В результате может суще- 
ственно снижаться частота вращения шпинделя, при которой наблюдается эффект параметриче- 
ского самовозбуждения. 
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Аналогичный эффект наблюдается и в том случае, когда в динамической системе резания 
рассматривается влияние двух физических процессов, вызывающих изменения параметров. Как и 
ранее, рассмотрим пример (рис. 3). В этом примере в отличие от (5) матрица жёсткости равна 


[1 + $1 (6. + 0,, ) + 5т(%,Ё + 3, ) [в эт( +0, ) + 5т(%,Ё + 3, > ) 
СЕС . (6) 
[4 эт( +8, >) + (0, + 3: > ) [1 + $ (6. + 0, >) + $т(0Ё + 3, > ) 


@) 
Здесь удобно ввести в рассмотрение коэффициент К, = Я изменение которого определяет со- 
1 


отношение частот параметрического возбуждения. Приведённые на иллюстрации данные получе- 
ны при р, =0,5. На приведённых иллюстрациях не показана высокочастотная область, как 


наименее интересная для рассматриваемых прикладных задач. 


П = 6, /6, П= 0, /@›% 
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Рис. 3. Изменение области параметрического возбуждения в плоскости (п, и.) при варьировании 


© 
частотного отношения К, =-—2:а— К, =0,2; 6— К, =0,9 
0), 


Здесь также обнаруживается существенное снижение нижней частоты вращения шпинде- 

ля, при которой наблюдается параметрическое самовозбуждение системы, по мере сближения 
частот параметрического возбуждения. В реальной динамической системе резания, имеющей 
большее число степеней свободы и условий, вызывающих периодические изменения параметров, 
картина параметрического самовозбуждения существенно расширяется. 
Заключение. Параметрические эффекты в динамической системе резания играют заметную 
роль в формировании области устойчивости траекторий формообразующих движений. Если сле- 
довать традиционному представлению о повышении устойчивости, то по мере увеличения скоро- 
сти резания, связанной с возрастанием частоты вращения шпинделя, область устойчивости долж- 
на расширяться. На практике это не соответствует данным экспериментов. По мере увеличения 
частоты вращения шпинделя, начиная с некоторого её критического значения, система вновь те- 
ряет устойчивость. Это обусловлено параметрическим самовозбуждением. Поэтому в зависимости 
от состояния станка, вызывающего кинематические возмущения траекторий, и его конструктив- 
ных особенностей, по мере увеличения скорости резания существует ограниченный частотный 
диапазон вращения шпинделя, при котором в низкочастотной области теряется устойчивость 
процесса резания за счёт влияния, прежде всего, запаздывающих аргументов. В высокочастотной 
области устойчивость вновь теряется за счёт параметрического самовозбуждения системы. 
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РАВАМЕТВТС $ЕТЕ-ЕХСТТАТТОМ ОЕ СОТТ1МС ОУМАМТС 5У$ТЕМ” 


У. Е. ФаКоуогонту, РВат ТВи Ничопд 
(Роп Зае Тесптка! Упмег$Ку) 


Тре ргоМет оп а БИ=тд {те {оораЁй депегайоп ге/анме {ю {те игогкреесе {акта о ассоипЕ {пе рагатеЁ ис 5е!!- 
ехсйайоп /5 сопу!аегеа. Еасюгз саизтд {Те сусс рагатеег иапаНоп$ п {те сита аупат зузет аге апа/угед. 
ТИезе !асогх аге геаЁе юЮ !1е зрайа! ато{тору о! е/а5Нс ргорег@ез т те игогкр/есе зибзу&ет, а/оигапсе оп 
игогкр/есе гойавоп рептеег иагайоп$, Ктетайс ребиграйоп$ гот {те аме вагаиаге оЁ пе тасВте ехесийуе 
е/етепЕ, апа {0 ше репосбГс ргосеззез т пе си тд гопе. Тйе депетс аупатс тосе! о! {те репобГс зу&ет 15 
О/Меп. ТВе еТасЕ оЁ {пе репобса/у уагутд рагатеег5 оп {пе ра ау 15 ше!еа. ВисКИпа сойарзе Геаигез 
5рес!с юг сита адупат/ зубетб аге 05с/оед. 1 15 або пои {таЁ {Пгоидй {те рагатеётс ас@опз, ий! т- 
сгеазе т сита рева ипаег пе рта! гои/ дгоийИ, Ёпеге а№иауз ехибЁ5 сиИса! тедиепсу аЁ ириср те зубет 
105е5 5 заРИКу. 

Кеуигогаб&: 5а Шу, сус/с рагатеег иапайоп5, рагатейтс зе/Р‘ехсйайоп. 


* ТНе гееагсй {5 аопе мт {Пе Нате оЁ {Не паерепдепе В&р. 
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